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■ We propose a unified magnetic phase diagram of cuprate superconductors. A 

. new feature of this phase diagram is a broad intermediate doping region of quantum- 
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critical, z = I, behavior, characterized by temperature independent T1T/T2G and 
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^ . linear TiT, where the spin waves are not completely absorbed by the electron-hole 
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continuum. The spin gap in the moderately doped materials is related to the sup- 
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CJ . pression of the low-energy spectral weight in the quantum disordered, z=l, regime. 



The crossover to the z = 2 regime, where T1T/T2Q — const, occurs only in the fully 
doped materials. 
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1. Introduction 

Recent measurements of the spin-echo decay rate, I/T2G, for a number of cuprate 

oxides, taken together with earlier measurements of spin-lattice relaxation rate 1/Ti, provide 
considerable insight into their low frequency spin dynamics. In this communication we show 
how these measurements may be combined with straightforward scaling arguments to obtain 
a unified magnetic phase diagram for the Y- and La-based systems. 

In the presence of strong antiferromagnetic correlations at a wavevector Q, the main 
contribution to both Tj"^ and for copper comes from small q = q — Q, so that one may 
write 1^: 

-— ~/c/qlim , /c?qx(q,0) 

III J i^^O UJ J-2G I--' 

where ^(q, ^) is the electronic spin susceptibility near Q. On making use of quite straight- 
forward scaling arguments, when applicable, one may substitute x(q, cj) = ^'^^^x{ciC,,uj/uj) 
into Eq.(0) {—1] is the scaling dimension of the real space spin correlator) and obtain: 

^-Tr'^u^-^ ^r^e-", ^-ioj. (2) 

Jl ^20 J-2G 

where ^ is the correlation length and uj an appropriate energy scale. 

2. Applicable Scaling Regimes 

We consider first a clean quantum antiferromagnetic insulator (referred to as "insulator" 
hereafter), described by the 5=1/2 Heisenberg model with the exchange coupling J. Be- 
cause the spin stiffness, ps — 0.18 J 0, is small compared to J, the quantum- critical (QC) 
scaling regime [0], where the only energy scale is set by temperature, ujr^T, exists over a 
substantial temperature range 2ps^T J The dynamical exponent, z, which relates 

the characteristic energy and length scales according to ci}~^~^, is z = l as a consequence of 
Lorentz invariance at T = 0, reflected in the linear dispersion relation of the spin waves. In 
this case, T1T/T2G ~ cj^ ~ ~ const. Since ^ ~ Cu~^^'^ ~ y-i/z ^ 7^-1^ further obtains 



(1) 
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l/Ti~T^^ "^r^T^ and 1/T2g~^^ ^ r-^T^ because the critical exponent i] is negligible, 
1/Ti~ const, while 1/T2g^T~\ 

A second regime of interest is the two-dimensional renormalized classical (RC) regime, 
T/v ^ T ^ 2ps, which is characterized by an exponential increase of the correlation length 
and relaxation rates. In the dynamical scaling theory of Chakravarty, Halperin, and Nelson 
0, 1/Ti ~r3/2exp(27rp,/r) and I/T2G exp{2npjT). The prefactors arising from the 
log^ corrections lead to a power-law temperature dependence of the ratio T1T/T2G 
while z = l leads to the cancellation of the leading (exponential) terms. 

According to numerical calculations for the insulator in the 2D 5=1/2 Heisenberg 
model |[TO|JTT[| , as long as T< J, the damping, 7^, of the high energy (c<jg>c^~^) spin wave 



excitations, is small throughout the Brillouin zone; hence, for both RC and QC regimes 
those can be treated as good eigenstates of the model. The dynamical susceptibility can 
then be well approximated as: 

X{q,uj)=4>g[ ^— ■ — ^ ) , (3) 

except near the origin, where the dynamics is diffusive as a consequence of total spin con- 
servation, and near the Neel ordering vector, Q = (7r/a,7r/a), where it is relaxational; for 

^q^a"^, the expression (H) is valid in both QC and RC regimes, where (pg ~ 1/q and 
cijg^cg. According to Ref. 0], the one-magnon neutron scattering intensity in the insulator 
is indeed well described by Eq.(|^) with ■yg<t^u!g. 

It has been conjectured in Ref. that the small doping as well as randomness related 
to it are not likely to affect the universal scaling behavior at high temperatures. Quite 
generally, one would expect a departure from z = l behavior only when spin waves become 
overdamped by the electron- hole continuum |13[. Since this would require a substantial 
increase in spin wave damping, to '~^g>ujg, compared to its insulator value, 7q<^co'q, there 
may be an intermediate regime in which the spin waves are not yet absorbed by the electron- 
hole continuum, even if the damping is increased compared to the insulator. To the extent 
this occurs, the system can remain in the QC regime with z = 1 in a wide range of doping 



and temperatures due to the Lorentz invariant terms in the action. 

In the insulator, the zero temperature energy gap A = hc/^ for the spin-1 excitations in 
a quantum disordered (QD) regime is, again, related to the Lorentz invariance at T = [0. 
Hence, as long as the Lorentz invariant terms in the action are still important and the corre- 
lation length saturates, the low frequency {u < A) spectral weight could be suppressed even 
in a metal. We suggest that the spin gap phenomenon in YBa2Cu306.63 and YBa2Cu408, 
characterized by a sharp increase in TiT and decrease of the bulk susceptibility, is related 
to this suppression, and this phase corresponds to a QD, z = l, regime. 

At larger hole densities, the spin waves will be fully damped by the electron-hole con- 
tinuum. In this regime, the self-consistent renormalization (SCR) approach developed by 
Moriya et al. (see also Ref. [|13[) and the phenomeno logical theory of Millis et al. 



||15|| , for YBa2Cu307 may be expected to apply. The dynamical exponent z = 2, while 
?7 = 0; in contrast with the previous case, the mean field exponent 2; = 2 is not fixed 



by a symmetry, but rather follows from the scaling analysis of Ref. |14| , p!3| . One obtains 
TiT/T2^Q ~a}^^ ~ const, while T1T/T2G At high temperatures, the energy scale 

d) ~ T, so that ^ ~ o)""'^/^ ~ 

T-^h ^ T-^l^, in which case T1T/T2G ~ T^/^, and, separately, 
1/Ti ~T^^~ const, l/r2G~^ r^T^^/"^. We emphasize that 1/Ti ~ const at high tempera- 
tures is predicted for both z = l and z = 2 regimes, while predictions for I/T2G are different. 
Finally, at still larger hole densities, the short range AF correlations between spins will 
tend to disappear; in this limit, is independent of temperature, and one recovers the 

Korringa law, 1/Ti ~T, while I/T2G ~ const. This regime corresponds to a normal metal, 
in which any remaining antiferromagnetic correlations can be described by a temperature 
independent F"(p,p'). 

3. La2_2:Sra^^Cu04 

We first consider La2-xSr2;Cu04. The insulator, La2Cu04, is well described by the 2D 
Heisenberg model with the nearest-neighbor exchange coupling J ~ ISOO-ft', except near 
or below T/v ~ 300 K induced by weak interplanar coupling [|l^]. A nearly temperature 



independent 1/Ti is observed in the insulating La2Cu04 above 650K 0], as expected in the 
QC, -2 = 1, regime The absolute value of 1/Ti 2700 sec~^ at high temperatures 
is in very good agreement with both expansion ^ and finite cluster |Jll| calculations 
for the Heisenberg model. Further, the ratio T1T/T2G measured in the insulator is 
nearly temperature independent in the broad range 450K<T<900K (Fig.|^). This is again 
what one would expect in the QC, z = 1, regime; we note that this behavior holds even 
in the region below 650 K, where 1/Ti and T/T2G separately deviate from constant values, 
apparently because T1T/T2G is insensitive to the magnitude of x(q, cj). 

A finite cluster calculation in the 5 = 1/2 2D Heisenberg model by E. Gagliano, S. Bacci, 
and one of the authors (A.S.) (Ref. [jlT| and this work), with no adjustable parameters used 
(hyperfine and exchange couplings were determined from other experiments, see Ref. [0), 
indeed yields a nearly temperature independent T1T/T2G — 4.3 • lO^K^^ for T> J/2:i^750K 
(Fig.0), in excellent agreement with the experimental result, 4.5 ■ lO^K p|. The systematic 
error of the finite cluster calculation, arising from the periodic boundary conditions as well 
as spin diffusion contribution to 1/Ti, is estimated in Ref. ||TT| as 10 — 15%. 

An especially striking feature of the 1/Ti data is the nearly doping independent 
absolute value of 1/Ti in the high temperature limit (Fig.^). This result shows that not 
only universal, but also nonuniversal scaling constants are not strongly renormalized in the 
doping range x = — 0.15, i.e. up to the optimal concentration for the superconductivity. 
Since it would be rather unrealistic to assume that exactly the same value can be obtained 
in different pictures of magnetism for low- and high-doped La-based materials, we suggest 
that the high temperature magnetic behavior over this entire doping range has the same 
physical origin as that found for the insulating state, which implies z = 1. As the doping 
increases, Imai, Slichter, and collaborators find that the range of temperatures where TiT 
is linear in temperature stretches towards lower temperatures, from 650 K for a; = down to 
125 K for x = 0.15 (see Fig.^. This is the behavior expected if doping leads to a decrease 
in ps, thus extending the QC region We thus conclude that Lai.85Sro.i5Cu04 is in the 
quantum-critical, z = l, regime for T^125K. 
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As it is evident from Fig.^ (inset), for 125K, TiT begins to depart from its linear in T 
behavior, exhibiting an upturn for Tf5i60 — 70K. We attribute this effect to the suppression 
of the low frequency spectral weight in the quantum disordered, z = l, regime [|I^. At lower 
doping values, experiment shows that the low temperature phase is not a superconductor but 
rather a spin glass, in agreement with the scaling analysis [Q, while for the lowest doping, 
the low temperature phase is the antiferromagnetic Neel state. 
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4. YBazCugOe 63 and YBa2Cu408 



YBa2Cu306.63 and YBa2Cu408 have quite similar properties. The product TiT measured 
in YBa2Cu306.63 [0)0]; is hnear in temperature for 160K ^ T ^ 300K, while it exhibits 
similar behavior in YBa2Cu408 for 170K^T^800K. Since a linear TiT is predicted in both 
quantum-critical (-2 = 1) and overdamped {z = 2) regimes at high temperatures, to distinguish 
between these regimes, we turn to the I/T2G data on YBa2Cu306.63 0, and plot T1T/T2G and 
T1T/T2G as a function of temperature (Fig.|D. In the range 200K<T<300K, T1T/T2G is 
nearly constant, while T1T/T2Q varies significantly, in agreement with the prediction for z = l. 
Were this material in the z = 2 regime, T1T/T2G would increase as the temperature increases, 
while T1T/T2Q would be constant. We thus conclude that above 200K, YBa2Cu306.63, and 
the closely related YBa2Cu408, are in the QC, z = l, regime. The increase in damping in the 
doped case, which enhances 1/Ti with respect to I/T2G, may explain the smaller (compared 
to the insulating La2Cu04, Fig.g) saturation value of T1T/T2G in those compounds. Our 
scenario may seem to contradict the Raman studies in YBa2Cu306+x, because two-magnon 
Raman scattering is not observed for doping above 0^ 4 |]TU|. However, this contradiction is 
illusory, since the decrease of intensity in the two-magnon Raman scattering is primarily due 
to the loss of the charge-transfer states rather than any change in the short range magnetic 
correlations ||T9|j20[ . 

For temperatures below 150 K, 1/Ti sharply drops down as the temperature decreases, 
while I/T2G 0,0 saturates. As was the case for Lai.85Sro.i5Cu04, we argue this suppres- 
sion of the low frequency spectral weight {spin gap) reflects a crossover to the quantum 
disordered, z = 1, regime. Since in this regime the magnitude of the gap is inversely pro- 
portional to the correlation length, which is smaller in YBa2Cu306.63 and YBa2Cu408 than 
in Lai.85Sro.i5Cu04, the onset temperature of the quantum disordered (spin gap) regime is 
larger and the crossover to it is more pronounced in the former two materials. Thus, unlike 
the scenario proposed by Millis and Monien [|T^], we argue that for all three materials the 
physical origin of the spin gap is the same. 
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5. YBaaCugOr 

For nearly stoichiometric YBa2Cu306.9, T1T/T2G increases as the temperature in- 
creases, while T1T/T2Q is nearly constant (Fig.^, in agreement with the scaling prediction 
for the overdamped, z = 2, regime The departure from the Korringa law 1/Ti ~ T 



and large copper-to-oxygen ratio of the relaxation rates shows that the antiferromagnetic 
enhancement is still quite substantial. In the overdamped regime, the spin wave branch 
is either destroyed, or due to the small correlation length has appreciable spectral weight 
only for energies much larger than the maximal temperature at which this compound is 
chemically stable. Therefore, for experimentally accessible temperatures no departure from 
the the overdamped regime is observed. 

6. Conclusion 

We have shown that the nuclear relaxation data in a broad range of doping levels, which 
includes both metallic and insulating materials, possesses universal features characteristic 
of the quantum critical regime of a clean antiferromagnetic insulator with the dynamical 
exponent z = l. This universality suggests that well beyond the metal-insulator transition, 
the spin excitation spectrum of the insulator is not yet destroyed by the electron-hole back- 
ground, in which case one expects a two-component dynamics for a broad range of doping 
levels. While a two-component dynamics may arise in a one-component as well as a two- 
component microscopic model, we call attention to an explicit example which leads to this 
kind of dynamics directly, namely, a model of spins and fermions with both spin-spin (J) 
and weak spin-fermion (J) exchange interaction [|21[. The robustness of the spin waves may 
be related to either their weak coupling to quasiparticles, or to quasiparticle Fermi surfaces 
which are not spanned by the antiferromagnetic ordering vector Q. 

On the basis of the above analysis, we suggest the unified magnetic phase diagram for 
the cuprate superconductors shown in Fig.|T]; the proposed boundary between the QC and 
QD regimes is determined from the nuclear relaxation data shown on Fig.^ We propose 
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that as the hole doping increases, the transition from the insulating to the overdoped regime 
occurs in two stages. First, the system becomes metallic; the damping of spin waves in- 
creases somewhat compared to its value in the insulator, but since the spin waves are not 
destroyed by the electron-hole background, the dynamical exponent is z = 1 and the quan- 
tum critical regime persists over a wide range of temperatures and doping levels. Then, 
at substantially higher doping, the dynamical exponent crosses over to z = 2. We further 
argue that the spin gap phenomenon observed in the underdoped materials reflects the same 
physics as the formation of the gap for spin excitations in the the quantum disordered, z = l, 
phase of an insulator. This scenario suggests that in compounds where the spin gap is ob- 
served, the temperature dependent bulk susceptibility should exhibit a downturn near the 
crossover from the QC to QD regime, while T1T/T2G should be temperature independent 
at higher temperatures. We show, in a subsequent communication [^, that our scenario 



leads in a natural way to a unified description of the results of nuclear relaxation, magnetic 
susceptibility, and neutron scattering experiments. 
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FIGURES 



FIG. 1. The proposed magnetic phase diagram for the cuprate superconductors above lOOK. 
FIG. 2. Experimental data on TiT and TiT/Tsg: La2Cu04 |||; 



V Lai.85Sro.i5Cu04 §; A Lai.85Sro.i5Cu04 |2|; • YBa2Cu408 [H; o YBasCusOg.es Hi; ■ 
YBa2Cu30e.9 Also shown (O) are the results of numerical calculation of 1/Ti and I/T2G 
( |11| and this work) for the insulator. The arrows indicate our proposed values for the crossover 



temperature from the QC to QD regimes; the inset makes clear our choice for Lai.85Sro.i5Cu04. 

FIG. 3. Experimental data on T1T/T2G and TiT/T^^: o YBa2Cu306.63 |1J|; 
■ YBasCugOe.g 
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